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Abstract 

We study time-dependent perturbations to a family of five-dimensional black hole spacetimes 
constructed as a holographic model of the QCD phase diagram. We use the results to calculate two 
transport coefficients, the bulk viscosity and conductivity, as well as the associated baryon diffusion 
constant, throughout the phase diagram. Near the critical point in the T-[i plane, the transport 
coefficients remain finite, although their derivatives diverge, and the diffusion goes to zero. This 
provides further evidence that large- N c gauge theories suppress convective transport. We also find 
a divergence in the low-temperature bulk viscosity, outside the region expected to match QCD, and 
compare the results to the transport behavior of known R-charged black holes. 
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1 Introduction and Summary 

The phase diagram of quantum chromodynamics (QCD) is of great interest. At vanishing chemical 
potential for baryon number, the theory is expected to evolve from a gas of hadrons at low tem- 
perature to a plasma of liberated quarks and gluons at high temperature. Were chiral symmetry 
exact, this evolution would progress through a true phase transition; given the breaking of chiral 
symmetry by quark masses, however, it is predicted by lattice gauge theory instead to be a smooth 
but rapid crossover (for recent lattice results, see [H [2] [3]). The results of heavy ion experiments 
at the Relativistic Heavy Ion Collider (RHIC) and now at the Large Hadron Collider (LHC) are 
consistent with such a picture (see for example [1]). However, near the critical temperature, a gas 
of quarks and gluons is not observed: instead, the physics is well-described by the hydrodynamics 
of a near-perfect fluid, with no apparent quasiparticle description. Such a fluid is a good candidate 
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to be modeled by a gravity dual via the AdS/CFT correspondence (for reviews, see [5j El [7]), and 
indeed the observed very small ratio of shear viscosity to entropy density [8] is a generic feature of 
such models [91 [10] . 

For QCD at nonzero chemical potential, it is expected that the crossover sharpens into a line 
of true first-order phase transitions, ending on a critical point. The temperature-baryon chemical 
potential (T-fj,) phase diagram is illustrated in figure [TJ for reviews see |11[ \12\ I13j . Based on 
dimensionality and symmetry, this critical point is expected to be in the universality class of the 
three-dimensional Ising model, like the liquid-gas transition of an ordinary fluid. Although it has 
not been accessed experimentally so far, the critical point could be probed in future heavy ion 
experiments at RHIC, LHC or the proposed fixed-target Compressed Baryonic Matter (CBM) 
project at the Facility for Antiproton and Ion Research (FAIR) 115] . It is difficult, however, 
for lattice gauge theory techniques to be used at nonvanishing chemical potential, due to the sign 
problem of the fermion determinant. While various lattice techniques such as reweighting |16[ I17j 
and imaginary chemical potential [18] have been employed, and other methods such as Nambu- 
Jona-Lasinio models have been used |19L 1201 [21] , another theoretical approach would be welcome, 
and the AdS/CFT correspondence offers such an approach. 

In |22j . we constructed a holographic model of the QCD critical point. A bottom-up method was 
used: the gravity dual was constructed using the minimal set of fields that captured the essential 
dynamics. These fields were a metric tensor characterizing the geometry, a scalar field to encode 
the running of the gauge coupling, and a gauge field dual to the U(l) current for baryon number, 
to provide the chemical potential; note that from the AdS/CFT point of view, going to nonzero 
chemical potential involves merely adding one new field. Points in the phase diagram correspond 
to black hole solutions with varying Hawking temperature and electric charge. 

The resulting five- dimensional gravitational theory contained some freedom: the scalar potential 
and the gauge kinetic function, both of which were arbitrary functions of the scalar field. To 
capture the behavior of QCD, these functions were chosen to make the zero-chemical potential 
thermodynamics, including quark number susceptibility, match the crossover behavior predicted 
by lattice QCD, building on work in (231 [23]; a related body of work includes [25, 26J and the 
review |27j . While the potentials used in our work do not come from any known string theory 
construction, they were chosen as combinations of exponentials to match the qualitative features 
of five- dimensional supergravity. 

Having thus completely constrained the theory by matching to lattice QCD on the T-axis, 
ensembles of charged black holes were obtained numerically to fill out the phase diagram in the 
T-fj, plane. As anticipated, the crossover indeed sharpens into a line of first-order phase transitions 
at finite fj,; black holes corresponding to both phases as well as the intermediate thermodynamically 
unstable configuration were found. Moreover, the static critical phenomena of the critical point 
were studied, and the critical exponents were found to be consistent with the Ising mean field 
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Figure 1: The expected phase diagram of QCD. The line ending in a star is the first-order chiral 
transition and its critical endpoint. Below is the nuclear matter transition. At lower right are color 
superconducting phases, color-flavor locked and otherwise. 

values (a, /3,7, 5) = (0, 1/2, 1, 3). The location of the critical point in this model was predicted to 
be (T c ,/i c ) = (143 MeV, 783 MeV). Thus a realistic, consistent holographic realization of the low-^ 
QCD phase diagram emerged Q 

It is interesting to move beyond static phenomena, and to study dynamics. Time-dependent 
critical behavior includes transport properties, relaxation times and the response to time-dependent 
perturbations, all necessary for the understanding of a heavy ion collision evolving near the critical 
point. There are three relevant transport coefficients: the shear and bulk viscosities associated to 
the traceless and trace perturbations of the energy- momentum tensor, fj and C, respectively]^] and 
the conductivity A of the U{\) baryon density. The ratio of shear viscosity to entropy density is 
known to be universal fj/s = l/4ir in all two-derivative gravity models [9] [10], but A and ( are 
expected to depend nontrivially on the location in the phase diagram. Moreover, the dynamic 
critical exponent z controls the phenomenon of critical slowing down, where the equilibration time 
r grows with the correlation length ^ asr~f; this behavior is expected to determine how close 
to criticality a heavy ion collision can approach in the limited time available before freeze-out and 
hadronization, as ^ < (time) 1/7 2 [29] 130] . 

Static critical phenomena are sorted into universality classes by properties such as their dimen- 
sionality and symmetry of the order parameter. This concept was extended to dynamic critical 
phenomena by Hohenberg and Halperin [31] , who classified various universal "models" based on 
their conserved quantities, which manifest as hydrodynamic modes, and the Poisson brackets be- 
tween them and with the order parameter; dynamic universality classes are then determined by 

1 Expected color superconducting phases at high chemical potential cannot appear in this simple model, and are 
an interesting task for the future; for progress see for example 28 . 

2 In keeping with the literature, we denote the shear viscosity fj to avoid confusion with the critical exponent r\. 
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these models as well as the static class. Son and Stephanov argued that in QCD, only one com- 
bination of the baryon density and the chiral condensate survives as a hydrodynamic mode, and 
thus having this single conserved mode as well as conserved energy-momentum, the theory should 
fit into dynamic model H [32]. 

The response of a fluid to charge inhomogeneities is controlled by the diffusion constant D, as- 
sociated to the dispersion relation oj = —iDk 2 . The diffusion constant is related to the conductivity 
A and the charge susceptibility \ by 

D = ± (1) 

X 

and at a critical point, D tends to zero while x diverges. The behavior of A near the critical point, 
however, depends on the dynamic universality class. Hydrodynamic models neglecting nonlinear 
interactions predict that A remains finite at criticality; this is characteristic of Hohenberg and 
Halperin's model B, where the single hydrodynamic mode is taken to be a conserved density. The 
only kind of conduction possible in such a model is via diffusion. Model B predicts the dynamic 
critical exponent 

z = A-rj, (2) 

where r/ is the usual static exponent giving the anomalous dimension of the density two-point 
function. In model H, on the other hand, the inclusion of energy and momentum as hydrodynamic 
modes makes convective conductivity possible, and this new channel naively dominates. In model 
H, the conductivity A and shear viscosity fj diverge at the critical point, 

\~\T-T C \*>, fj~\T-T c \ x i, (3) 

with exponents x\ and x^ which are related to the static exponent r\ by 

x\ + x v = 4- d-ri, (4) 

where d is the spatial dimensionality. In the 3D Ising model rj and x^ are close to zero, giving x\ 
close to one [32J. The critical exponent z when convection dominates takes the value 

z = 4 - 7] - x x , (5) 

and is thus moved from z ~ 4 in model B to z ~ 3 in model H. Note that even in models where A 
diverges, \ always diverges faster, so D still goes to zero [H] . 

The expected behavior of bulk viscosity near the liquid-gas critical point was investigated by 
Onuki [33], who predicted a divergence also depending on z, 

\T-T c \~ zu+a (6) 
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where v > and a are static exponents; for several other predictions for bulk viscosity at a critical 
point, see for example [53] . 

Work has been done applying the AdS/CFT correspondence to dynamic critical phenomena in 
various models. Bulk viscosity for the models we consider has been computed at vanishing chemical 
potential in [35, 24J. Maeda, Natsuume and Okamura calculated the conductivity in the "one-charge 
black hole" model dual to M = 4 Super- Yang Mills with a chemical potential — hereafter the one- 
charge N = 4 black hole — and found it to stay finite at the critical point, identifying the result with 
model B despite the presence of conserved energy and momentum [36] . The natural speculation 
is that the nonlinear interactions responsible for the convective component in the conductivity in 
model H are suppressed by the large number of colors N c , and in [57] Natsuume and Okamura argued 
that large--/V c counting indeed enhances the diffusive over the convective conductivity, reducing 
model H to an effective model B. Meanwhile, the bulk viscosity was studied in a mass deformation 
of the one-charge M = 4 black hole by Buchel [38] and in the N = 2* model by Buchel and 
Pagnutti [39] [3D]. In both cases the bulk viscosity was also found to be finite at the critical point, 
in contradiction with expectations from ([6]). It is natural to speculate that again the large- N c limit 
suppresses the divergence in the transport coefficient. The dynamic critical exponent z has also 
been studied directly, yielding the mean-field model B value z = 4 for the deformed one-charge 
black hole [38] and z = for N = 2* [UJ. 

Dynamic universality classes depend not just on the conserved quantities and their Poisson 
brackets, but also on the static universality class. The work described in the previous paragraph 
shares the list of conserved quantities with QCD, but differs in the static critical exponents. It is 
natural to wish to extend the study of holographic dynamic critical phenomena to a model that 
shares the static universality class with QCD as well. In this paper, we begin the investigation 
of the dynamic critical phenomena in the QCD-like holographic critical point of [22J, by studying 
finite-cu fluctuations around those black hole backgrounds. With these fluctuations, we are able to 
calculate the transport coefficients C and A, and the associated diffusion D. 

In keeping with the other AdS/CFT cases, we find the transport coefficients remain finite at 
the critical point, and the diffusion constant goes to zero. Thus these models share the apparent 
suppression of the convective contribution to transport by the large- N c limit [37J, and behave 
effectively as model B. 

Furthermore, the behavior of A and £ as functions of T and p near the critical point is quite 
similar to that of the entropy and baryon densities s and p: all are smooth approaching the critical 
point in the T-p plane along the axis defined by the first-order line, but develop infinite slope when 
approaching off-axis. Similar behavior arises in the one-charge M = 4 black hole, where the critical 
exponent controlling all these divergent slopes is the same, 

A - A c ~ s - s c ~ p - p c ~ \T - Tel 1 / 5 , (7) 
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with 5 = 2; the corresponding exponent in the QCD-like case, 5 ~ 3, is consistent with our results 
for A — A c and C ~ Cc- This suggests that in the vicinity of the critical point, the deviations of the 
conductivity and bulk viscosity from their critical values can be thought of as depending smoothly 
on the deviations of the densities from theirs. 

From these results we can estimate the dynamic critical exponent z assuming a mean-field value 
of 77; a precise determination of z and r] requires finite-/c fluctuations, which we leave for future 
work. Given the model B behavior and the mean field exponents, using the mean field value r/ = 
in ([2]), we can estimate 

z«4. (8) 

Another feature of the behavior of the transport coefficients in the phase diagram is worth point- 
ing out. In [22j and [23], the potential and gauge kinetic function in the Lagrangian were only 
constrained to lattice QCD results over certain ranges of temperature, corresponding to certain 
values of the scalar field. In principle, one could imagine modifying these functions such that the 
matched region is not affected, while results change elsewhere; thus the application of these models 
to QCD should only be trusted in a certain band of temperatures. This is not unreasonable, since 
the gravity dual picture is only expected to apply in a certain range near the critical temperature 
where there is no quasiparticle description. 

However, phenomena outside the region designed to match QCD may be interesting in their 
own right. We find that as the temperature decreases, both the bulk viscosity and conductivity 
begin to rise. In fact, the bulk viscosity has a divergence at a temperature around half the crossover 
temperature on the T-axis, outside the region matched to QCD, and this divergence extends out into 
the plane. Unlike the phenomena at the critical point, this divergence in the transport coefficient is 
not associated with any feature in the thermodynamics of p or s. From the gravity point of view, it 
can be understood as the place where the fluctuation solutions develop a node, but the field theory 
interpretation is unclear. Whether the conductivity also has a pole at a still lower temperature is 
outside the range accessible to our numerical solutions. We note that a similar divergence appears 
for the conductivity of the thermodynamically unstable branch of the one-charge black hole, as the 
"superstar" limit [32] is approached. 

The summary of the remainder of the paper is as follows. In section[2]we discuss the Lagrangian, 
equations of motion and gauge symmetries of our class of models, and construct the gauge-invariant 
fluctuations around the black hole backgrounds. The fluctuation equations and the Kubo formulae 
for extracting the transport coefficients are presented in section [3} In section |4j we review the one- 
charge M = 4 black hole solutions and solve for the conductivity both analytically and numerically, 
checking our method and pointing out several features that will have analogs in the QCD-like black 
holes; the analytic calculation was previously worked out in [361 143] . In section [5j we present the 
calculations of the transport coefficients for the QCD-like black holes. We conclude in section [6] 
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Results for the conductivity and diffusion of another model, the two-charge TV = 4 black hole, are 
given in the appendix. 



2 Black hole backgrounds 

In this section we review the five-dimensional gravity ansatz that applies both to the QCD-like 
holographic critical models of [22] and the one-charge N = 4 black hole, present the resulting 
equations of motion and symmetries, and construct the gauge-invariant fluctuations around these 
backgrounds with nonzero frequency. 

2.1 Lagrangian and ansatz 

We consider five-dimensional gravitational theories containing a metric g^ v , a vector field and 
a scalar (j), with Lagrangian 



2k 2 



(9) 



The relevant solutions to these theories are asymptotically AdS black three-brane solutions with 
radial profiles for the electric potential and scalar field, 

ds 2 = e 2A{r)^_ h{r)dt 2 + df2) + ^^ dr 2 

Hn (io) 

A^dx^ = $(r)dt, (j) = (j)(r), 
and the equations of motion resulting from this ansatz are 

A" - A'B' + V 2 = 
h" + (4A' - B')ti - e- 2A f{(j))& 2 = 

„ + („_, + *),_£<£_o. 

where 

V«&, r) = V(<P) - I e -^-2B /(0)$ /2 ^ (12) 
along with the zero-energy constraint 

h(2AA' 2 - (j)' 2 ) + GA'ti + 2e 2B V{(j)) + e - 2A /(</>)$' 2 = . (13) 
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Asymptotically AdS solutions may be written in coordinates that as r — > oo approach 

ds 2 -> e 2r/L (-dt 2 + dx 2 ) + dr 2 , 
<S>(r) = <S> {0) + $ {2) e- 2r / L + ... , (14) 
Hr) = ^(4-A)e (A " 4)r/L + . • • + <P {A) e- Ar/L + .... 

Black brane backgrounds have a horizon r = th defined by the largest solution to the vanishing of 
the horizon function h{rn) = 0. The temperature T and entropy density s can be calculated as 

T = L h >( rH ) e MrH)-B(r H ) g = 2 _l e ZA { r H ) (15) 

while the chemical potential [i and £7(1) density p can be read off from the near-boundary expansion 
of <&(r), 

, = %i, ,= -*f. (16) 

2.2 Gauge symmetries 

The gauge symmetries of the theory are U(l) gauge transformations of the vector field given by 
X(x), and general coordinate transformations given by e^(x), with the general transformation of 
the fields 

5 = d^X + e v d v A ll + {d^e u )A u , 
54> = e»d^ . 



To stay consistent with the ansatz (10), one may still make transformations with a general e r (r) as 
well as e°(t), e l (x) and X(t) obeying 

dt€°(t) = const , diCj(x) + djCi(x) = const x 6{j , dtX(t) = const . (18) 

One may use e r (r) to choose the metric function B(r) arbitrarily. In the gauge B{r) = 0, which 
we use for the QCD-like holographic critical black holes, the remaining coordinate transformations 



(18) may be used to set any two of the zero point of A(r) and the overall scales of h(r) and $(r): 



(r,x\t) = ( -/r,f3x\-t ] -»■ h = \h, A = A + logf3, $ = -$ , (19) 
V 7 / 7 2 7 

while the gauge transformation A can fix the zero point of $(r). The remaining transformations are 
constant shifts of all variables (of which shifts of t and x % are trivial, while shifts of r will change 
the form of the various functions), and the manifest SO(3) symmetry acting on x % . 
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2.3 Fluctuations and gauge invariant quantities 

We are interested in studying small fluctuations around the background ( 10 ) of the form 



ds 2 



5° + Re{e 2A ^e-^% u (r)} dx^dx" , 

L J (20) 

A^dx!* = $(r) dt + Re{e- Wt a lt (r)}dx li , <f> = <f>(r) + Re{e~ lujt <j>(r)} , 



where g® u is the background metric, and we have assumed plane- wave dependence in the time 
direction. This ansatz preserves the SO (3) rotating the spatial directions together. From here on, 
we will omit the instruction to take the real part and work directly with complexified perturbations, 
on the understanding that appropriate superpositions of complexified perturbations will be real. 

Certain combinations of the fluctuations will be gauge degrees of freedom. These correspond 
to gauge transformations at linear order in the small parameter associated to the fluctuations; 
gauge transformations at zeroth order act on the background. Gauge transformations preserving 



the ansatz ( 20 ) for the time dependence of the fluctuations also take a plane wave form, 

e"(t, x, r) = e~ iwt e(r) , A(i, x, r) = e' iljt r]{r) . (21) 

We can then calculate the transformation of the fluctuations huv, (p under the gauge transfor- 
mations. Here we record only the results for modes with no index in the r-direction; all independent 
fluctuating modes come from this set of fields. For the metric we have, 

6h tt = -{ti + 2A'h)C + Muhf , 

5hu = -iuj5ij(, j , (22) 
5hij = 2A!CSij , 



for the gauge field, 



and for the scalar, 



$a t = -iujri + f r & - icL><3?£° , 

(23) 

Sat = , 



We would like to construct gauge-invariant combinations corresponding to the fluctuating modes of 
the theory. We expect a total of nine such fluctuations: five modes from the graviton, three from the 
gauge field and one from the scalar. We immediately see eight such gauge invariant combinations: 
the three <Xj, and the five traceless components of hij. The trace of the spatial part of the graviton 
transforms as 

1 



5 ( 3 (h xx + h yy + h zz ))= 2A'C , (25) 
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and so we can construct the ninth gauge invariant mode involving the scalar fluctuation, 

S = 4> - ^^(h xx + h yy + h zz ) . (26) 

Thus the nontrivial profile for the background scalar field (ft = <j)(r) couples the scalar fluctuations 
to fluctuations of the graviton trace; in the AdS limit, this reduces to just the scalar fluctuation 
S — > (ft and the graviton trace becomes a pure gauge mode. For a nonzero scalar profile, this is the 
AdS/CFT encoding of the running coupling that breaks scale invariance producing a nonzero trace 
of the energy- momentum tensor. 

While S is a natural definition for the fluctuation since it has a smooth AdS limit and shares 
its possible asymptotic behaviors with (ft, it will be useful for us to consider the rescaled mode 

l 2A< (27) 

= -(h XX + hyy + k ZZ ) ~ — (ft . 

This normalization is the natural one for thinking of the mode as the graviton trace; for (ft 1 ^ one 
may pick a gauge where (ft = 0, where T~L is the trace precisely. 

We see the fields organize themselves under SO (3) as a singlet (the scalar S or %), a triplet 
(the gauge field a^) and a quintuplet (the traceless graviton hij), and since they are in different 
representations they cannot mix at the linearized level. Thus we expect each of these modes to 
satisfy a decoupled fluctuation equation. Note this would not hold for a nonzero spatial momentum; 
instead, a smaller SO(2) symmetry group would organize the perturbations. 

3 Fluctuation equations and transport coefficients 

The fluctuation equations of the nine modes 7~L, ai and hij can be obtained by linearizing the 



Einstein, Maxwell and Klein-Gordon equations around the background (10). In general, the various 
equations that result will also contain the auxiliary quantities a r , h^r- By taking linear combinations 
of the Einstein and Maxwell equations, these can be eliminated, giving us nine differential equations 
in nine unknowns. 

The coefficients of the fluctuation equations include the background fields A, B, h, $ and (ft and 
their derivatives. These fields obey the background equations of motion, which act as constraints 



on these coefficients. To exhaust the constraints, we use the various background equations (11) to 
eliminate the second derivatives A", <3?", h" and (ft"; in addition, the zero-energy constraint (13) 
can be used to eliminate the potential V((ft), although derivatives of the potential will still appear. 
In this way we achieve a unique presentation for each coefficient with no hidden constraints. 

Each fluctuation equation can be used to calculate a transport coefficient via a Kubo formula. 
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These are associated to the imaginary part of the corresponding Green's function, which comes 
from the "conserved flux" T for each equation. For a differential equation of the form 



y"(r) + p(r)y'{r) + q(r)y(r) = , (28) 

Abel's identity guarantees that for two solutions yi(r), y2 (r), the quantity 

T = exp ^ j p(r)dr) W( yi ,y 2 ), (29) 

is independent of r, where W (2/1,2/2) = 2/12/2 ~~ 2/22/'i is the Wronskian of the two solutions. As long 
as p(r), q(r) are real, for some complex solution y(r) its conjugate y*{r) is also a solution. Taking 
2/1 = 2/ and 2/2 = 2/*! our expression for the conserved flux becomes 



T = exp y J p(r)drj Im (y*y) , (30) 

up to an overall factor. 

3.1 Near-horizon behavior and boundary conditions 

Noting that h has a zero at the horizon, h = /i'(rjj)(r — m) + . . ., while A and its derivatives are 
regular there, all the fluctuation equations considered in this paper have the same limit near the 
horizon, 

1 2 2B(r H )-2A(r H ) 

X" + — —X> + V2 X = , (31) 

where X is any of Z, a, S or T~L (see the following subsections), and we have assumed that the 
gauge degree of freedom B is chosen so it and its derivatives do not diverge at the horizon. Near 
the horizon, one may expand 

X(r) = (r - r H ) a (x + x x {r - r H ) + . . .) , (32) 

where we must choose the exponent a to be 

e B(r H )-A(r H ) 

a = ±iuj — u? — ^ — • ( 33 ) 

h'(r H ) 

In solving the fluctuation equations we must impose two boundary conditions. The first is the 
requirement of infalling boundary conditions at the horizon: this corresponds to solving for retarded 
Green's functions, in accord with the usual prescription for calculating transport coefficients. This 



simply involves imposing that only the negative sign solution in (33) contributes. 



The second boundary condition must be imposed not at the horizon, but at the boundary. Each 
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mode has two solutions near the boundary, one falling off more quickly and corresponding to a VEV 
deformation of the dual field theory, and the other falling off more slowly and corresponding to 
turning on a source. Our prescription will be to normalize the coefficient of the "source" mode to 
unity. 



Assuming we have imposed infalling boundary conditions, a part of the conserved flux (30) 
simplifies, 

hlm(X*X') = -uj\x Q \ 2 e B{rH) - Ai - rH) , (34) 

where we used the fact that a is pure imaginary. One can evaluate the conserved flux at any r, but 
since one boundary condition is at the horizon and the other at the boundary, one must in general 
solve for the mode everywhere in order to do so. As we shall see, in some cases it is possible to do 
this analytically in the cj — > limit; in general, we shall solve the fluctuation equations numerically. 

3.2 Traceless graviton and shear viscosity 

This case is well-known, and the outcome, fj/s = 1/4-7T, is guaranteed by the general arguments of 
[HUH]]. However, it is useful to go through an explicit calculation as a warmup for the more difficult 
conductivity and bulk viscosity cases, to be treated next. 

Letting Z be any of the traceless hij, the corresponding fluctuation equation is the massless 
scalar equation, 

/ U\ 2B-2A 

Z " + ( 4A' — B' + --) Z' + —r^- u 2 Z = 0. (35) 



h) h 2 
The behavior of Z near the boundary is 

Z = Z {0) + ... + Z {4) e- 4r / L + ... , (36) 

and so the proper boundary condition is simply to take ZtQ\ = 1. The associated conserved flux 
from Abel's identity is 

T z = he 4A - B lm(Z*Z'), (37) 
and the Kubo formula for the corresponding transport coefficient, the shear viscosity, is 

fj = -— lirn -T z . (38) 



The simplicity of the equation (35) in the uj — > limit allows one to solve for the shear viscosity 
analytically; we take our discussion from [Mj. In the us — > limit the term with no derivatives 
vanishes, and we are left with 

d r (log Z') = -d r (4A -B + logh), (39) 
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which has the solution 



poo -4A+B 

Z = a + b I dr . (40) 



The second term is technically not allowed at strict u = due to a logarithmic divergence; it may 
be kept for very small cj, but for us it is enough to note that matching to the near- horizon expansion 

Z(r) « z (r - r H ) a = z (l + alog(r - r H ) + . • .) , (41) 
we have Zq = clq; however as r — > oo we see ao = -Z(o)- Thus the boundary condition Zm\ = 1 



imposes Zq = 1 near the horizon. Using (34) and (37), one can then evaluate r/ at the horizon, 



fj = ^ MrH) , (42) 



which using ( 15 ) implies for the shear viscosity to entropy density ratio the familiar universal result, 

? = i-. (43) 

S 47T 

The other transport coefficients we consider will not be universal in this way; from a practical 
standpoint, imposing the boundary condition at infinity will not impose a universal constraint on 
the near-horizon behavior. A similar technique to the one reviewed here will, however, be useful in 
calculating the conductivity of the one-charge N = 4 black hole of the next section. 



3.3 Gauge field, conductivity and diffusion 

Because of SO(3) symmetry, we can treat each of the fluctuations separately: they cannot mix 
at linear order. For notational simplicity, we will drop the index i altogether and use a to denote 
one of the components of a^: for example, a = a% to compute the conductivity in the x 1 direction. 
The fluctuation equation is 

, + (W + fm ^ . _ „ , ,44, 

and the associated conserved flux is 



Ta = hf(<P)e 2A - B Im(a*a'). (45) 
Near the boundary the gauge field fluctuation behaves as 

a = a {0) +a (2) e- 2r/L + ... , (46) 
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so the appropriate boundary condition is simply = 1. The Kubo formula for this mode 
determines the (zero-frequency) conductivity A, 



A 



— 2 hm -F a 



2k 2 w->o oj 

This quantity is related to the frequency-dependent complex conductivity, defined as 

i L 0(2) 



(47) 



a 



oj a 



(o) 



by 



A = lim Re<r(u;) 



(48) 



(49) 



which follows from (45) and (48) using the boundary condition. The diffusion constant D is related 

X = D X , (50) 
dp" 



to the conductivity by 
where the susceptibility is 



X 



dp, 



leading to the the Nernst-Einstein relation, 



D\ = lim Recrfa;) . 



(51) 



(52) 



Unlike the shear viscosity, the conductivity does not take a universal value; notice that the equation 



(44) still has a zero-derivative term at u = 0. For the example of the next section, however, we will 



be able to solve the w — > limit and find an analytic solution for the conductivity. 

3.4 Scalar fluctuation and bulk viscosity 

Finally the scalar equation is 

S" + (V -B' + jjs'+ {^—^u 1 + E(r)) 5 = 0, 
where S(r) is determined by the background: 

n -2A 



(53) 



18fh 2 A' 2 



18^'Y V 2 - e 2A fh 2 ct>' A + GfhA'tf 



2e 2A+2B V 



+ e 2A ti(j)' + f'& 2 +3/A 



/2 B „2Ai2i)2 



/iV 2 + 3h& 2 f" - 6e 2A+2B hV" 



(54) 



The associated conserved flux is 



T s = he iA ~ B Im (S*S r ) , 



(55) 
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and one may in principle calculate the bulk viscosity £ from this flux using the formula 

C = -Ai im ij- 5j (56) 

with S suitably normalized. Since a gauge may be chosen where S = <f>, it is evident that 5 has 
the same asymptotic behavior as the scalar <fi itself, 

S(r) = S (4 . A) e^ A -^ L + ... + S (A) e- rA / L + .... (57) 



Thus to calculate (55), one must properly normalize the coefficient 5(4_a) °f the- dying exponential. 
This is not as straightforward numerically as normalizing an asymptotic constant. 

Since the bulk viscosity is a quantity extracted from the trace of the energy-momentum tensor, 
and the energy-momentum tensor couples to the graviton, it is the mode T~L that we more naturally 



wish to canonically normalize. The scaling of % depends on <p{r) via (27). If the background scalar 
profile includes the "source" term 0(4_a)> we have A' /(f)' ~ e r ( 4 ~ A )/ L and we end up with the % 
scaling 

H(r) = n {0) + ... + %_ 2A) e r ( 4 - 2A )/ L + . . . , (58) 

and the suitable boundary condition to calculate the bulk viscosity is %(o) = 1> an easier prescription 
to implement^] 

The fluctuation equation for % simplifies most when we use the background equations of motion 
to eliminate V and (ft' 2 in favor of </>" and A": 

H" + (aA! -B' + j + 2 -f- 2 -§-)^+ (^- 2 + ^(r)) H = 0, (59) 

with 

h' / A" 4>"\ . e- M ,„ .,2 



E " = J U-W) + -W (3A "'- ' <60) 

One can check that the boundary behavior as described above does indeed result from the asymp- 



totic expansion of (59), as it must. The coefficient of the first-derivative term is again explicitly 



integrable, and it is then easy to use Abel's identity to show that the conserved flux is 

P 4A-B y/2 

?U = - J? Im(7W), (61) 



where we chose the normalization so that this conserved flux coincides exactly with J 7 ^ (55) when 



3 The asymptotic scaling of H is different when the background scalar profile only contains the VEV term 4>(a) 
but no source. However, in this case the conformal invariance of the dual field theory is not explicitly broken, and 
the bulk viscosity is identically zero. 
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one substitutes (27): 

^=?s(s^-&). (62) 
Thus the bulk viscosity may most easily be calculated as 

C = -7T^ lim n -F-H, (63) 

with the boundary condition = 1- 
3.5 Large-A^ counting 

The various Kubo formulae are all proportional to the five-dimensional gravitational constant 1/k 2 , 
since they all come ultimately from evaluating the five-dimensional action. The gravitational con- 
stant can be expressed in terms of the AdS radius L and the underlying string coupling g s and 
Regge slope a' when an explicit string theory construction is specified. For AdS§ x S 5 , the result 
is 5 

= 647r 4 5s V 2 ' (64) 

while the AdS radius is 

L 4 = 4irg s N c a' 2 , (65) 

with N c the number of colors in the dual gauge theory (the five-form flux from the gravity point 
of view). Field theory quantities should not involve a' , which drops out of the combination 

T 3 N 2 

±^ = (66) 
k 2 4vr 2 ' 

and all the Kubo formulae indeed result in this combination on dimensional grounds. As a result, 
all three transport coefficients go like N 2 in the large- N c limit defined by the dual gauge theory. 
For AdS/CFT models not based on a known string theory construction, the precise coefficient 



in (66) is not determined. However, the N 2 dependence is expected to remain the same for any 
gravity dual of a four-dimensional gauge theory. For the QCD-like holographic critical models of 
section [5j we will use the formula (66) with the understanding that the overall normalization is not 



truly determined, but the N 2 factor is expected to be robust. Note too that the factor (66) cancels 
out of C/s and D, and so will only be relevant for us in X/T. 



4 One-charge J\f = 4 black hole 

Before discussing our primary interest, the numerical QCD-like critical black holes, we first consider 
a family of black hole backgrounds where the thermodynamics are known analytically, which we 
call the one-charge M = 4 black hole |44l 1431 IHH W7\ H"51 |4*U] . These geometries are known solutions 
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of string theory, coming most simply from a truncation of the maximally supersymmetric gauged 
supergravity in five dimensions, which is in turn a truncation of the dimensional reduction of type 
IIB string theory on AdS 5 x S 5 . The dual field theory configurations should be thought of as states 
in M = 4 super- Yang-Mills theory with a temperature and a chemical potential for a £7(1) subgroup 
of the SO(6) R-symmetry. In the appendix, we discuss a second family, the two-charge M = 4 
black hole, which does not display finite-temperature critical phenomena. 

Since the temperature and chemical potential are the only massive parameters, dimensionless 
quantities can only depend on their ratio; for this reason the phase diagram is not truly two- 
dimensional, but is more properly thought of as depending on this single ratio. The phase diagram 
has the form of a semi-infinite line, ending on a critical point. 

Because the conformal invariance of the theory is not broken explicitly, the bulk viscosity is 
identically zero. The conductivity and associated diffusion have been calculated by [36] H3]. We 
reproduce the calculation for two reasons. First, since it is analytically solvable, it allows us to 
check our numerical methods against a known analytic solution; and second, it exhibits several 
features that will persist to the QCD-like black holes. 

In this section, in keeping with the literature, we use fi to denote a parameter in the solutions, 
and thus use Q for the chemical potential. Also, it is inconvenient to set B(r) = for these 
solutions, so for this section (and the appendix) only we employ a radial coordinate where AdS 



space is related to that given in (14) as exp(r t h er e/^) = r hem/L 



4.1 Thermodynamics and phase diagram of the one-charge M = 4 black hole 

The one-charge N = 4 black hole is a solution to the Lagrangian Q with the potential and gauge 
kinetic function, 

V W) = ~J2 ( 8e ^ +4e"^) , f((l>) = e" 2 V^. (67) 
The solution takes the form 



r 1 / Q 2 \ r 1 / O 2 

A(r) = log - + - log 1 + ^r , B(r) = -log---log 1 + ^ 



L 6 b V r2 J ' L 3 

<■(>•) = 1- , «r) = -,/Ilog(l + 5T), .(0- ^ Q ^ 



r 2 (r 2 + Q 2 )' Yy ' V3 B V r 2 /' r 2 + Q 2 r\ + Q 2 



(68) 



and is characterized by the charge parameter Q, mass parameter [i and the asymptotic AdS scale 
L. The horizon is at 

th = J I (VQ 4 + 4A^-Q 2 ) , ( 69 ) 
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and in what follows we will generally trade the parameter fx for r#. The temperature and chemical 
potential can be expressed as 



Q + 2r| n= _ Qr H 



2nL 2 JQ 2 + r 2 H L 2 J Q 2 + r 2 H 



(70) 



For fixed Q and L, \x is only required to be nonnegative; the limit fx — > corresponds to rjj — > 
but the temperature approaches 

T(^0)^^. (71) 

This value is properly thought of as the limiting temperature of a sequence of honest black holes 
with nonzero values of fx. The solution with strictly fx = is not a black hole, as the horizon 
function is trivial h(r) = 1 and no horizon exists. The fx = solution is supersymmetric and has 
been dubbed the "superstar" by Myers and Tafjord [12]. 

Choosing Q to be positive along with rjj, T and Vt are in general multivalued as functions of 
these parameters. The transformation 



r H ^ ^ 5, Q^r H V g , (72) 

2 \ Q 2 + r 2 H \ Q 2 + r 2 H 



maps pairs (rn,Q) to physically distinct pairs with the same (T, Q). The transformation takes 

Q 2 , 2r| 



2^ Q2 

which defines a fixed locus, 



(73) 



Q = V2r H -»• vrT = \/2ft . (74) 

Since dimensionless quantities like s/T 3 and x/T 2 will be functions only of the dimensionless ratio 
O/T, we should properly think of this system as a one-dimensional phase diagram with this locus 
interpreted as a fixed point 0,/T = 7r/v2, rather than a fixed line. The models of the next section 
add another massive parameter, with explicitly breaks conformal invariance and leads to a true 
two-dimensional phase diagram. 
We notice that in general, 

which establishes the minimum value of T for a given f2 lies at the fixed point, 

T > — . (76) 

7T 
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Thus for vrT < there are no corresponding black holes, at the fixed point ttT = t/2Q there is 
one, and for ttT > \/20 there are two. The limit Q — > produces two configurations dual under the 



transformation (72): the superstar with r# = 0, Q 7^ 0, and an uncharged black hole with Q = 0, 

Black holes on the two branches with the same (T, O) will not have the same values of other 
thermodynamic quantities, and so are physically distinct. The entropy and charge density in terms 
of Q and th are 



s 



which in terms of T and Q can be written 



X ' 1 *- (3tt T 2/) V ± V) , P = T V? , (78) 



167T 167T 2 



where the top sign corresponds to the branch with > Q, the bottom sign to the branch with 

s/2th < Q, and where we have defined 

y = Vk 2 - 2& , (79) 



with = n/T. Note that the positivity of the quantity inside the square root follows from (76). 
The quantities s and p are the T and Q derivatives of the pressure, 

p = if? (" 4 + 5vr ^ 2 " T ± ^) ' (80) 

We note that the free energy density is / = —p, so for fixed (T, Q) the V^fH > Q branch has a 
lower free energy and hence is thermodynamically preferred in the canonical ensemble. 
The U(l) susceptibility for the two branches is 

* = f? (5^-30^^), (81) 

which diverges along the fixed line for y = 0. The heat capacity and off-diagonal susceptibility 
diverge as well. Thus the fixed point is a second-order phase transition. The determinant of the 
matrix of susceptibilities S is 

„ 3iV r 4 T 4 / 4 9 ~ 9 ~, o 46vr 3 ^ 2 - llvr^ 4 \ , x 

det S = — 4vr 4 - 22ir 2 n 2 + ft 4 ± 4ir 3 y ± , (82 

64^^ V J 

and while this is always positive for the thermodynamically preferred branch \p2rn > Q, it is always 
negative for the other; thus the non-preferred branch is thermodynamically unstable as well. 
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The exponents for this critical point have been calculated. The fundamental relation is the 
approach of the density to the critical density, which for fixed T takes the form: 



p-p c ~ \n-n c \ 1/2 = \n-n c \ 1/s , (83) 

defining the critical exponent 

5 = 2. (84) 

An analogous relation holds for the approach to T c with fixed f2, as well as for the entropy. The 



divergence of the conductivity x ( as wen as the heat capacity Cq) follows from (83), 

x~ |fi-o c |- 1/2 = \n-n c \~ e , (85) 

where the exponent]^] e = 1 — 1/5 = 1/2. 

4.2 Conductivity and diffusion of the one-charge J\f = 4 black hole 

We may now turn to solving for the conductivity. First we consider the gauge field fluctuation 



equation (44). This can be solved exactly at u = 0, allowing an analytic determination of the 
conductivity and the associated diffusion constant using similar techniques to those used for the 
shear viscosity. We find the uj = solutions 

Q 2 + 2r 2 

a(r) = Cl (Q2 + 2)(Q2 + r2) + C2 a 2 (r) , (86) 
where a>2{r) is a more complicated expression including log(r — rn) whose explicit expression is 



unenlightening. In order to determine C\ and C2 , one may match the zero- frequency solution ( 86 ) 
to the near-horizon solution 

a(r) = (r - r H ) a , (87) 



where the exponent a from (33) corresponding to infalling boundary conditions is 



iuL 2 ^Q* + r 2 H 
a = ~ 2Q 2 + Ar 2 H 



4 In the literature for critical phenomena, the well-known exponent labels a, f3 and 7 are often associated to 
behavior approaching the critical point along the axis denned by the first-order line. Since there is no first-order 
transition in this model, we avoid these exponent names, e encodes the approach of \ t° the critical point off the 
first-order axis. 
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Then looking at the near-boundary behavior, we can determine the series (46) and solve for the 



frequency-dependent conductivity (48). For small uj we find the result 



*Q 2 , L(Q2 + 2r ^2 
a ~ o — 27 + / + ° w i (89) 



which in terms of T and £1 can be written 



a N: 



'c 

T 4?r 2 



^ (2vr 2 - n 2 T 2vry) + + 0( W ) 

2a; V / 7r ± y 



(90) 



where = fi/T and u = oj/T. The 1/w pole in Im<r indicates by the Kramers- Kronig relations a 
delta- function contribution to the real part, Re a oc 5(to), not visible in our calculation. This infinite 
contribution is characteristic of a translationally- invariant charged system |50| . Setting aside the 
delta-function, we have the zero-frequency conductivity |36j 

4(vr±y) 

We also calculated this conductivity using the numerical methods that will be employed in the next 
section, and find excellent agreement. 

Since y — >• at the critical point, the conductivity A approaches the constant 

A c = thh . (92) 

4-7T 

The conductivity not diverging is a sign of the large- N c domination of diffusive conduction (charac- 
teristic of model B) over convective conduction (model H) [36, 37J. However, the slope of A diverges 
as the critical point is approached. One has 

A - A c ~ (T - T c ) 1 /2 . (93) 

We note that the exponent is 1/5; this transport coefficient has the same critical behavior as the 
thermodynamic quantities s and p. We will find an analogous phenomenon for the QCD-like black 
holes. 

The diffusion D = X/x is given by 

D= A 4 1 (94) 

Since the conductivity is constant at the critical point, while the susceptibility diverges, the diffusion 
goes to zero. On the stable branch at large temperature, D asymptotes to its AdS value D — > 1 /2irT. 
Because x on the unstable branch is negative close to the critical point and passes through zero, 
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Figure 2: The conductivity over temperature and diffusion times temperature for the one-charge 
black hole with $7 = 1, for the stable branch only. X/T approaches a finite value with an infinite 
slope at the critical point, while DT goes to zero. 
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Figure 3: The conductivity over temperature and diffusion times temperature for the one-charge 
black hole with Q = 1, for both branches; the unstable branch is the thick red line. The approach 
to infinity at large T is the superstar divergence. 



the diffusion on the unstable branch is negative and diverges before returning from positive infinity. 

Away from the critical point, an entirely distinct phenomenon occurs in the superstar limit. 
Since as y(p, — )■ 0) — )■ ir, the conductivity ratio X/T approaches zero at the charged black hole 
(the far end of the stable branch) and approaches infinity in the superstar limit (the far end of the 
unstable branch). The rate of approach to infinity in the superstar limit becomes 



A 



1 



(95) 



We will find a similar divergence in the bulk viscosity of the holographic critical black holes of the 
next section, away from the critical point and the region expected to be related to QCD. 
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5 QCD-like holographic critical black holes 



We now turn to our primary interest, studying the transport coefficients of the class of solutions 
designed to emulate the thermodynamic and phase structure of QCD |22| . These black holes are 
solutions to the Lagrangian ansatz ([9]) with the potential 

= -12cosh70 + ^ 2 with = Q 606 and b = 2m7 (g6) 



and the gauge kinetic function, 

, , sech [| w 
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The black hole solutions are of the form (10). For convenience the gauge B(r) = was chosen, and 



the residual coordinate transformations ( 19 ) and the freedom to shift r were used to set 

A(r H ) = 0, *(rjr) = 0, ti(r H ) = l/L, r H = 0. (98) 
The solutions were characterized by two initial conditions at the horizon, <pQ = </>(r#) and $i = 



$'(rjf); the ensemble of black holes was generated by numerically integrating the equations (11) 
from the horizon to the boundary for a large set of distinct (<^o,$i), and the thermodynamic 
quantities of temperature T, entropy density s, chemical potential fi and baryon density p were 
obtained for each. 



The functions (96) and (97) were chosen so that the thermodynamics of the fi = black holes 
reproduce lattice results for the equation of state [23] and the quark susceptibility [22], building 
into the model the rapid crossover that connects the hadron phase to the quark-gluon phase in 
physical QCD with massive quarks; see figure [4| Normalizations for the scales of T and fj, were also 
derived from matching to the lattice results, allowing the thermodynamic quantities to be expressed 



in MeV. The potential (96) leads to the scalar mass, 

m\L 2 « -0.293 , (99) 

with implies a dual operator that is barely on the relevant side of marginal, 

A^« 3.93, (100) 

designed to emulate the slow running of the QCD coupling constant somewhat above the confine- 
ment scale. The profile of the scalar in the background includes a "source" term that explicitly 
breaks conformal invar iance. 

At nonzero p for this ensemble of black holes, the crossover sharpens into a line of first-order 
phase transitions ending at a critical point, as is expected for QCD. The region of the first-order 
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Figure 4: The normalized entropy s/T 3 and quark susceptibility \2 = X2/T 2 at p, = 0, computed 
on the lattice and fit by black holes in the gravity theory defined by our choices of V((f>) and f(4>) 
(equations (96) and (97)). Lattice data is taken from p]. 



line is characterized by the existence of two distinct phases, realized as distinct black hole solu- 
tions at the same point in the phase diagram, as well as a third black hole corresponding to the 
thermodynamically unstable state lying in between the stable phases in the Maxwell construction. 
The first-order line ends on the critical point. For the static critical exponents defined as 

/ ds\ 

C p = T ( 7^ ) ~ \T — T c \ 01 , along first order axis 



dT J 



p 



Ap ~ (T c — T)' , along first order line 

X = ( -J- ) ~ \T — T c |~ 7 , along first order axis 



(101) 



dn J 



T 



p- p c ~ \n- fi c \ 1/S , forr = T c , 

this critical point was found numerically to have 

a = 0, /3ss 0.482, 7 w 0.942 , 5 rj 3.035 , (102) 

with accuracies consistent with the Ising mean field values (a, /3,7, 5) = (0, 1/2, 1, 3). Matching the 
temperature and chemical potential scales to that of the lattice data, the location of the critical 
point was found to be 

T c = 143 MeV p c = 783 MeV . (103) 

While these results generate a large phase diagram, the applicability to QCD is expected to be 
limited to a band surrounding the crossover and extending out into the plane, for two reasons. 



First of all, the potential (96) and gauge kinetic function (97) were only matched to lattice data 



over a finite range. In principle, one can imagine varying the functions so as to keep matching 
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the data in the range required, while changing the results elsewhere. The potential V(4>) was only 
matched up to about <po « 7.5, corresponding to T m i n rj 135 MeV at \i = 0; the crossover is at 
T c ~ 175 MeV on the T-axis. Moreover, AdS/CFT is only expected to provide a good description 
of QCD in the region where no quasiparticle description is available; both above and below the 
crossover such a description is possible, with hadrons at lower T and quarks and gluons at higher T. 
Thus we expect our geometries to be reasonable potential models for QCD only near the crossover; 
and we hypothesize that this applicability extends along with the crossover out into the T-fi plane 
to the critical point. 



5.1 Boundary conditions 

To calculate the correct transport coefficients (or thermodynamics for that matter) for these black 
holes, a technical detail must be taken into account. The numerical integration beginning at the 
horizon produces solutions that asymptote to AdS space, but in general in coordinates that do not 



match standard AdS coordinates. In the gauge B(r) = 0, one usually uses the AdS coordinates ( 14 ) 



with A(r) —7- r/L and h(r) — > 1 + 0{e~ Ar / L ); however with the coordinate choices (98) used for the 
numerical solution, the horizon function in general approaches a non-unit value h — > , and the 
warp factor goes like A(r) — > r/{\Jh^L) + A^ 1 ' . To calculate thermodynamic and hydrodynamic 
quantities using the usual formulae, it is necessary to pass to a new set of coordinates that takes 
the usual AdS form, with unit normalization for the horizon function and A ~ r/L; this can be 



done by making a different rescaling (19). 



Moreover, the scalar field generically approaches <f) —> <f)A^~ uA ^ r \ with v = 4 — ~ 0.07. It 
is useful to rescale the leading perturbation to the scalar to a standard magnitude <j>A = 1, which 
can be achieved by adding a constant to r. Fixing the value of 4>a is useful because turning on 
this perturbation corresponds to deforming the ultraviolet theory with an almost-marginal relevant 
operator, playing the role of the running coupling, and hence changing the theory. We want every 
point on the phase diagram to correspond to the same theory, not different theories, and hence 
we make (f>A identical in all of them. Put another way, there are three massive parameters in the 
theory, T, fi and A, where A is the scale associated with the scalar. Normalizing (j)^ = 1 corresponds 
to measuring T and /i in units of A. The explicit breaking of conformal invariance coming from A 
permits a non-zero bulk viscosity to be present. 

The desired coordinates are achieved by the transformation 



i = $ V h$* t, ~x = <t>fx , T - = — T =- + At - log(<^ 



L ~ fZ^rr ° ~ 1 ' ( 104 ) 





which is a combination of a scale transformation associated to ^Y" and a time dilation by * / h^ r 
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as well as a shift of r, implying the relation of the functions, 

r,_, 1 



A(r) = A(r) - log^J ) h(r) 



1 



^farMO <*>(r) = -— =$(r). 



A V 



(105) 



Thermodynamic and hydrodynamic quantities are naturally calculated in terms of the tilded co- 
ordinates because of their standard near-AdS form. The black hole solutions exist in the untilded 
variables, but as we now show it is possible to compute certain ratios directly in the untilded co- 



ordinates. One may also use the thermodynamic formulae (15), (16) in the untilded variables to 



calculate "horizon" quantities Tjj, [in-, sh and pn- The true entropy and temperature are related 
to their "horizon" counterparts as 



1 



SH 



(106) 



where the horizon quantities are constants independent of the black hole considered, 



2tt 



4ttL 



SH 



(107) 



Transforming from one set of coordinates to the other, and using uj = <$~a \i h^ui, we have 5 



Thus, we can assemble the invariants 

Vh _ V 

SH S ' 



Xh 



\h_ 

S H 



,1/3 



Ch- 



Ch = C 

SH S 



and thus these ratios can be calculated directly in the untilded variables. 



(108) 



(109) 



With the gauge choices B{r) = and (98), the value of the exponent a characterizing the 



near-boundary behavior of the fluctuation equations is simply 



a = ±iLcu . 



;no) 



As usual, one should choose the minus sign to impose infalling boundary conditions at the hori- 
zon. The boundary conditions at the boundary are that a and % approach unity. As mentioned 
previously, if one employed S instead of H, it would be the coefficient of the e~ vr term that had 
to be properly normalized; this is numerically more difficult. The results outlined in the following 



5 Although under this rescaling one has S = a/cfy , we assume the gauge field fluctuation is asymptotically 
normalized to one in both coordinate systems before being plugged into the Kubo formulae. 
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Figure 5: The bulk viscosity over entropy density over the T-fjL plane for the QCD-like black hole 
solutions. 



subsections were obtained using the H equation, but the S method was checked as well, and while 
the S calculation developed numerical pathologies in certain (low-T) regions of the phase diagram, 
the two methods agreed well in the regions where both appeared reliable. 



5.2 Bulk viscosity for the QCD-like holographic critical black holes 

The computation of the bulk viscosity for the QCD-like black hole solutions is a straightforward 



application of the Kubo formula (63). To construct the mode solutions to the % fluctuation 



equation (59) we used a near-horizon series expansion to seed a numerical integration routine 
which propagates the mode solution out to close to the boundary. 



Because the flux (61) satisfies Abel's identity, it must be independent of the radial location at 
which it is evaluated; ascertaining that this holds is a check that the fluctuations obtained indeed 



satisfy (59). After verifying that this remains true in practice, we then evaluated the flux at a 
radial location where the numerical solutions are particularly robust. The black hole solutions 
are defined from r ~ 10 -6 near the horizon out to r ~ 10, and we used a value around rx ~ 10 . 



With the fluctuation in hand, we normalized it to one asymptotically and inserted it into (63) and 
used (106) and (109) to calculate the bulk viscosity over entropy density ratio C/ s 
6 Since fj/s is a constant, £/s is simply proportional to Q/f). 
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Figure 6: The bulk viscosity over entropy density at zero chemical potential and at [i = [i c for the 
QCD-like black hole solutions; the "bump" at left evolves into the peak with singular slope on one 
side as n is increased. 

The results of this calculation, when carried out for many points on the T-\x plane, are shown in 
figure [5] up to jjL = jjL c ; in all figures in this section we have used the normalization of |22| to express 
T and [i in MeV. Two features are immediately apparent. One is the propagation of a "bump" on 
the T-axis at the location of the crossover, out into the plane towards the critical point. As the 
bump moves into the plane, it becomes increasingly asymmetrical as the slope on the lower-T side 
becomes more vertical. Precisely at the critical point, this slope diverges. Constant-^ plots of £/s 
for both fj, = and [i = [i c are given in figure [6j The result at vanishing fi was obtained previously 
in [21]. For \i > /i c , the peak "tips over" and the plot of C,/s becomes multivalued. 

The first and more important conclusion from this critical behavior is that the bulk viscosity 
remains finite at the critical point; it does not diverge as is predicted in a number of models such 
as that of Onuki [33] , given in equation ^ ; see also [33] . 

Moreover, the way in which the quantity stays finite is interesting. Approaching the critical 
point along the "crest" of the peak turns out to be equivalent to approaching along the axis defined 
by the first-order line (but on the other side of the critical point); along this path, the value of (,/s 
scarcely changes at all. On the other hand, approaching the critical point from a direction other 
than the first-order axis, for example constant fi = yu c as depicted on the right of figure [6j Q/s 
develops a divergent slope; the same is true for constant T or any other direction of approach, save 
only the first-order axis, where no such divergent slope is seen. 

The bulk viscosity is not the only quantity with this property. It is common for thermodynamic 
functions to behave differently depending on whether the approach to the critical point is on 
the first-order axis or not. In [22] it was found that the thermodynamic densities s and p also 
behave smoothly approaching the critical point along the first-order axis, but have divergent slope 
approaching off-axis. Since the derivatives of these densities are the specific heat and susceptibilities, 
these properties are characterized by various static critical exponents: the smooth approach along 
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Figure 7: The bulk viscosity over entropy density at zero chemical potential down to lower tem- 
peratures, showing the divergence. 



the first-order axis corresponds to the vanishing a = 0, while the divergent slope approaching off 
the axis is encoded in the value of 5 « 3, giving an approach of the susceptibility off the first-order 
axis as 



X~ \fi- He 



with 



\T-T r 



off first order axis , 



2/3. 



Ill) 



(112) 



A natural guess is then that this behavior can be explained by the hypothesis that the bulk viscosity 
is a smooth function of the two densities: C, = p). Let us characterize the divergent slope of C, 
by an exponent 5^: 



C-Cc 



\T-T c \ 1/5 < 



;ii3) 



Then we would expect 6^ = 5. We analyzed the results of numerics and found an exponent 
consistent with this interpretation, although the strong sensitivity to the precise value of £ c , which 
is difficult to determine due to the infinite slope, gave an uncertainty in the exponent in the range of 
0.1 < l/#£ < 0.7. A similar argument can be carried through more precisely for the zero-frequency 



conductivity A of the one-charge M = 4 black hole. In this case, the result (93) shows that 5\ = 2, 
in agreement with 5 = 2 for the one-charge model. 

The second salient feature of the Q/s on the phase diagram is the strong rise at colder temper- 
atures and smaller chemical potential. This behavior is shown for zero chemical potential in figure 
[7| which shows the same thing as the left-hand-side of figure [6] but including lower values of T. The 
bulk viscosity actually diverges at a temperature T*, which on the T-axis is at T*(/i = 0) ~ 86 MeV, 
before becoming finite again at smaller temperatures; this divergence extends from the T-axis out 
along a curve into the finite chemical potential region of the phase diagram. The divergence is 
power law, with 



(T — T*) 



;ii4) 
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Notably, the thermodynamic densities s and p and their derivatives do nothing special at this locus. 

From the gravity point of view, this divergence can be understood as related to a formation of 
a node in the mode solution for %{r). For T > T*, the fluctuation %{r) never crosses zero as it 
asymptotes to a constant value at the boundary. For T < T* on the other hand, the solutions have a 
single zero before asymptoting to a constant. As T — > from below, this node moves out towards 
the boundary, and precisely at T = T* the solution goes to zero at infinity; that is, the leading 



constant in the asymptotic expansion of % (58) vanishes. Since the AdS/CFT prescription is 



to normalize this constant to unity, the scaling required to accomplish this normalization leads to 
the divergence. 

A reasonable expectation is that for T < T*, there is a perturbative instability in the black hole. 
Because the thermodynamics is perfectly well-defined and stable in the vicinity of T = T*, such an 
instability, if present, would be a decisive counter-example to the correlated stability conjecture of 

On the field theory side, however, it is unclear what this divergence means, and what the 
conjectured instability would be for T < T*. Importantly, this particular feature is not a prediction 
for QCD-like systems, since it lies well outside the range that was matched to lattice QCD data; 

rs 86 MeV is well below the lower bound of T m i n ~ 135 MeV where the matching to lattice 
QCD ended. (In fact, the value of the temperature where Q/ s begins to rise towards the divergence 
is close to T m ; n .) Instead one should view this phenomenon as an indication of what features are 
generally possible in gauge/gravity transport. 

Interestingly, this divergence is reminiscent of the superstar divergence in the conductivity for 



the one-charge N = 4 black hole (95). In both cases, the transport coefficient has a double pole in 



the thermodynamic variable as the divergence is approached. 

5.3 Conductivity and diffusion for the QCD-like holographic critical black holes 

Utilizing a numerical algorithm identical in spirit to that employed for the computation of the 



bulk viscosity, it was also possible to construct solutions to the fluctuation equation for a{r) (44) 



and subsequently evaluate the conductivity Kubo formula (47). The results of this calculation are 
displayed in figure [8j where the dimensionless ratio X/T is plotted over the T-p plane. Unlike the 
ratio £/s, this quantity grows in the larger-T region of the diagram. Note that as mentioned in 
section [3~5l the overall normalization for this ratio is not known due to the lack of a string theory 
embedding, but the N% dependence should be general, so we plot X/TN^ using (66) as a concrete 
choice. 

Once again, the most salient attribute is the propagation of the crossover feature out to the 
critical point, which is the sudden rise of X/T. Cuts at \i = and \x = \i c are displayed in figure |9j 
As with the bulk viscosity, the conductivity attains a finite value at the critical point; and also 
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Figure 8: The conductivity over temperature over the T-fj, plane for the QCD-like black hole 
solutions. 
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Figure 9: Conductivity over temperature at zero chemical potential and at \i = \i c for the QCD-like 
black hole solutions; the smooth rise on the left evolves to the singular jump on the right as fi is 
increased. 
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Figure 10: Baryon diffusion at zero chemical potential for the QCD-like black hole solutions. 

like the bulk viscosity, the slope of the rise increases as the critical point is approached, until it 
diverges, except along the first-order axis. Thus it is again natural to hypothesize that A = A(s,p), 
and that off the first-order axis the conductivity approaches its critical value with exponent l/5\ 



A -A, 



\T -T r \ l/5X 



;i!5) 



with 5\ = 5. 

As described in the introduction, the finiteness of the critical conductivity is reminiscent of 
dynamic critical model B rather than model H [31], matches the result for the one-charge J\f = 4 

), and is consistent with the interpretation of a large- N c 



models (see equation ( 92 ) and [361 H 
enhancement of diffusive conductivity over convective conductivity |37j . 

We may also consider the baryon diffusion Db = A/%. Since the conductivity diverges at 
the critical point, the finite critical conductivity implies vanishing baryon diffusion. In [22J, the 
behavior of the baryon susceptibility was studied as the second order point was approached along 



the first order line axis, and along directions off this axis. From (102) and (111), one has that on 



and off the axis the susceptibility diverges with critical exponents 7 ~ 1 and e ~ 2/3, respectively. 
Thus, we find the rates of approach to zero for the diffusion 



T IT 



\T-T r 



along first order axis , 



;ii6) 



and 



IT — T r 



, r _ r ,2/3 



off first order axis . 



;ii7) 



In figure 10, the conductivity divided by temperature is plotted for vanishing chemical potential. 



In this region of the phase diagram, the baryon susceptibility is particularly easy to compute, as it 
can be explicitly related to the background metric and scalar [22J. As a result, it is straightforward 
to extract the baryon diffusion along the T-axis. The diffusion also experiences a rise at low 
temperature, but we are unable to determine whether or not this rise culminates in a divergence 
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(as is the case for the bulk viscosity) as it manifests in a numerically unreliable region of the phase 
diagram. 

6 Conclusions 

The primary lesson so far of the study holographic critical phenomena has been the suppression of 
fluctuations and convective transport by large N c . The models considered here have a thermody- 
namic crossover designed to emulate QCD, and correspondingly a critical point arises in the T-fi 
phase diagram. However, rather than exact 3D Ising static critical exponents, we found mean field 
Ising critical exponents; and rather than divergences in the transport coefficients, we here find finite 
critical behavior. (In fact, one can view the famous constancy of the shear viscosity over entropy 
density, which holds throughout the phase diagram, as being the first example of this phenomenon, 
since mode-coupling theories predict a weak divergence.) These results can be explained by the 
hypothesis that large- N c suppresses quantum corrections as well as convective transport |37j . 

Physical QCD of course has N c = 3, but large- iV c counting is still a useful way of looking at 
the theory for a number of phenomena. Understanding better the extent to which the suppression 
phenomena that emerge from the AdS/CFT duals apply to real QCD is the primary outstanding 
question: in other words, could large- -/V c -suppression in real QCD significantly push the behavior 
of physical quantities toward mean field and model B expectations, except very close to the critical 
point? Understanding both the nonzero-momentum fluctuations and the finite- N c corrections to 
these models is essential for making further progress, and we hope to report on these issues in the 
future. 
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Appendix: Two-charge M = 4 black hole 

Along with the one-charge N = 4 black hole solution of section |4j another family of black holes 
coming from string theory is the so-called two-charge M = 4 solution, which has equal charges for 
two U{1) gauge fields inside the SO(6) R-symmetry. Keeping only the diagonal gauge field, this 



33 



geometry solves our ansatz with potential and gauge kinetic function, 



V(<P) = -^[Se^+Ae V^j , /(«£) = e V§* , (118) 

where the scalar potential matches that for the one-charge case, while the gauge kinetic function 
is slightly different. The solution takes the form 

A{r) = \og r - + \\og(l + ^ , jB (r) = -log^-^log(l + ^) , 



— _ _ (119) 



(r 2 + Q 2 ) 2 ' ^ w V 3 6 V y ' v ' r 2 + Q 2 r\ + Q 2 ' 

This solution again is characterized by a charge Q, a mass parameter ^ and the asymptotic AdS 
scale L. The solution has a horizon at 

r H = ^L y fji-Q 2 , (120) 
For fixed Q and L, the mass parameter /j, is bounded below by 

and at /it = ^ m i ra we have r# = T = 0; extending /x < fj, m in results in a naked singularity. 

It is more convenient to trade the parameter fi for rn, which runs over all nonnegative values. 
The temperature and chemical potential are 

giving a one-to-one relationship between the black hole parameters (r#, Q) and the thermodynamic 
parameters (T, fi); a single black hole exists at each point in the phase diagram, and hence a single 
phase. The entropy and charge density are 

K ^ 4 2 (123) 



k 2 L^^ H > 8vr 2 
which can be obtained as derivatives of the pressure 



P=^(2vr 2 r 2 + 2 ) 2 . (124) 
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The U(i) susceptibility is 



X 



dn) 



T 



8vr 



2 (2vr 2 T 2 + 3ft 2 ) 



(125) 



and it along with other derivatives of s and p such as as the specific heat are everywhere well- 
behaved; this example has no phase transitions. 

We can consider the conductivity and related diffusion coefficient for this model. Proceeding 
analogously to the one-charge case, we find the oj = solutions 



a(r) = d 



Q 2 + r< 



+ C 2 a 2 (r) 



(126) 



where, as for the one-charge case, a 2 (r) is a more complicated expression including log(r — r#) 
Matching to the near-horizon solution 



a(r) = (r — rn) a 



(127) 



with the exponent a imposing infalling boundary conditions given as (33) 

iujL 2 



a 



(128) 



we determine C\ and C 2 and solve for the conductivity (48). For small oj we find the result 



a(u)) 



+ 



JQ 2 

ojk 2 L ' 2k 2 (Q 2 + t 2 h ] 



+ 0(ui). 



(129) 



It is straightforward to convert this to dependence on the thermodynamic quantities T, 0, and one 
finds 



a N 2 in 2 



- + 



7T" 



T 4vr 2 \ 2u 2vr 2 + & 



+ 0(u) , 



(130) 



where again 0, = f2/T, Cj = uj/T. As in the one-charge case we find al/w pole in the imaginary 
part, indicating a delta function in the real part, and the remaining zero-frequency conductivity 



A 



N 2 ttT 



4(2tt 2 T 2 + ft 2 ) 4(2vr 2 + ft 2 ) 



(131) 



which using the susceptibility ( 125 ) gives the diffusion constant 

1 



D 



1 



2nT Vl + 2ft 2 /7r 2 + 3ft 4 /47r 4 



(132) 



We plot X/TN 2 and DT in figure 11 
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Figure 11: The conductivity over temperature and diffusion times temperature for the two-charge 
black hole with $7 = 1. 
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